Distribution amplitudes and decay constants for (tt, K, p, K*) mesons in light-front 

quark model 
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We present a calculation of the quark distribution amplitudes(DAs), the Gegenbauer moments, 
and decay constants for it, p, K and K* mesons using the light-front quark model. While the quark 
DA for 7r is somewhat broader than the asymptotic one, that for p meson is very close to the 
asymptotic one. The quark DAs for K and K* show asymmetric form due to the flavor SU(3)- 
symmetry breaking effect. The decay constants for the transversely polarized p and K* mesons(/J 
and ) as well as the longitudinally polarized ones(/ p and fx* ) are also obtained. Our averaged 
values for fy/fv, i-e. (/J//p)av = 0.78 and (/^» //K*)av = 0.84, are found to be consistent with 
other model predictions. Especially, our results for the decay constants are in a good agreement 
with the SU(6) symmetry relation, /J (x *) = (/vr(K) + f P (K*))/ 2 - 



I. INTRODUCTION 



Hadronic distribution amplitudes (DAs) are important 
ingredients in applying QCD to hard exclusive processes 
via the factorization theorem 0, B 13- They provide 
essential information on the nonperturbative structure of 
hadron describing the distribution of partons in terms of 
the longitudinal momentum fractions inside the hadron. 
Both the electromagnetic form factors at high Q 2 and 
the B-physics phenomenology are highly relevant to the 
detailed computation of hadronic DAs. 

During the past few decades, there have been many 
theoretical efforts to calculate the pion DA using non- 
perturbative methods such as the QCD sum rule [EBLIE 
US 11, lattice calculation [1, [l(| fill. fl^ . fl3ll , chiral quark 
model from the instanton vacuum 

[M UM, lli, Nambu- 
Jona-Lasinio(NJL) model [Vn. ITa|. and light-front quark 
model(LFQM) [l^, H3]. It is well known that the shape 
of the pion quark DA is very important in the predictions 
of pion electromagnetic form factor in both nonperturba- 
tive and perturbative momentum transfer regimes. The 
QCD sum-rule based analysis |g, 0] of the tt — 7 tran- 
sition form factor F 71 " 1 " 7 (Q 2 ) measured by the CLEO 
experiment (2lj has shown that neither double-humped 
DA for the pion predicted by Chernyak and Zhitnit- 
sky nor the asymptotic one are favored at the 2a 
level of accuracy. It is also interesting to note that 
the recent anti-de Sitter space geometry/ conformal field 
theory (AdS/CFT) prediction [22| for the meson DA is 
aVdS /cft {%) \J x(l — x), which would approach to the 
asymptotic form x(l — x) only in the limit of In Q 2 — ► 00. 
The shape of 0AdS/CFT( x ) increases the usual perturba- 
tive QCD (PQCD) predictions for the pion form factor 
and 77 — 7 transition form factor by 16/9 and 4/3, respec- 
tively. In our recent LFQM application to the PQCD 
analysis of the pion form factor p3| , we further found 
a correlation between the shape of quark DA and the 
amount of soft and hard contributions to the pion form 
factor. Similar to the previous findings from the Sudakov 
suppression of the soft contribution (or enhancement of 



the hard contribution) [2J, [2f| |2fJ l23| , our results indi- 
cated that the suppression of the endpoint region for the 
quark DA corresponds to the suppression(cnhanccment) 
of the soft (hard) contribution. 

Another important area that requires detailed study 
of meson DAs is the B-physics phenomenology under in- 
tense experimental investigation at BaBar and Belle ex- 
periments. The K, p and K* DAs have attracted at- 
tention rather recently [H, [2^, HH due to the deep 
relevance to the exclusive B-meson decays to (K, p, K*) 
mesons. In particular, the SU (3) flavor symmetry break- 
ing effect in the meson quark DA including strange quark 
is important for the predictions of exclusive i? Uj d iS -decays 
to light pseudoscalar and vector mesons in the context 
of CP-violation and Cabibbo-Kobayashi-Maskawa quark 
mixing matrix studies. The SU(3) breaking effect is real- 
ized in the difference between the longitunal momenta of 
the strange and nonstrange quark, (x s — x u ^) ^ 0, in the 
two particle Fock components of the meson. The similar 
effect was also found in our PQCD analysis [12] for the 
exclusive heavy meson pair production in e + e~ annihila- 
tions at a/s = 10.6 GeV. Not only the shape of the heavy 
meson quark DAs matters in the prediction of the cross 
section for the heavy meson pair productions, but also the 
cross section ratios for c(e + e~ — > DfD~)/<r(e + e~ — > 
D+D-) and cr(e+e- -> B^B° s )/a(e+e~ -> B+B~) de- 
viate from 1 appreciably due to the SU (3) symmetry 
breaking. 

A particularly convenient and intuitive framework 
in applying PQCD to exclusive processes is based 
upon the light-front (LF) Fock-state decomposition of 
hadronic state. In the LF framework, the valence quark 
DA is computed from the valence LF wave function 
ty n (xi, kj_i) of the hadron at equal LF time r = t + 
z/c which is the probability amplitude to find n con- 
stituents(quarks,antiquarks, and gluons) with LF mo- 
menta ki = (xi,h±i) in a hadron. Here, Xi and kj^ are 
the LF momentum fraction and the transverse momenta 
of the «th constituent in the n-particle Fock-state, respec- 
tively. If the factorization theorem in PQCD is applicable 
to exclusive processes, then the invariant amplitude M. 
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for exclusive process factorizes into the convolution of the 
process-independent valence quark DA <fi(x, p) with the 
process-dependent hard scattering amplitude Th Q, i.e. 



[dxi] [dyi]<f>(xi,fi)TH(xi,yi,fx)if>(yi,fi), (1) 



M 



where [dxi] = <5(1— E^^x^n?, dx^ and n is the number 
of quarks in the valence Fock state. Here, p denotes the 
separation scale between perturbative and nonperturba- 
tive regime. Since the collinear divergences are summed 
in (f>(xi,fj,), the hard scattering amplitude Th can be 
systematically computed as a perturbative expansion in 
a s (p) . To implement the factorization theorem given by 
Eq. |T]) at high momentum transfer, the hadronic wave 
function plays an important role linking between long 
distance nonperturbative QCD encoded in DA and short 
distance PQCD encoded in Th- 

The quark DA of a meson, c/)(x,p), is the probability of 
finding collinear quarks up to the scale p in the L z = 0(s- 
wave) projection of the meson wave function defined by 



|kJ<M 



[d 2 k ±l }^{ Xi ,k A 



where 



[d 2 k A 



2(2tt) 3 5 



3=1 



m 



d 2 k^ 
-2(2tt) 3 



(2) 



(3) 



Simple relativistic quark-model based on the LF frame- 
work has been studied for various mesons (l9l. I20I l33l |34| . 
Although the proof of duality between the LFQM and 
the first principle QCD is not yet available, we have at- 
tempted to fill the gap between the model wave function 
and the QCD-motivated effective Hamiltonian [35|, l36l |. 
The essential feature of our LFQM (3f| [36| is to treat the 
gaussian radial wave function as a trial function for the 
variational principle to the QCD-motivated Hamiltonian. 
We saturate the Fock state expansion by the constituent 



quark and antiquark, i.e. H„ 



Hq + Vi n t, where the 



interaction potential Vint consists of confining and hy- 
perfine interaction terms. From the variational princi- 
ple minimizing the central Hamiltonian with respect to 
the gaussian parameter, we can find the optimum values 
of our model parameters and predict the mass spectra 
for the low-lying ground state pseudoscalar and vector 
mesons [H, [36[. We applied our LFQM for various ex- 
clusive processes such as the electromagnetic form factors 
of 7r, K and p [35l Wf\ mesons and semileptonic and rare 
B decays to tt and K [36|, [38[ . Our results for the above 
exclusive processes were in a good agreement with the 
available data as well as other theoretical model predic- 
tions. 

The purpose of this work is to calculate the quark 
DAs, the Gegenbauer moments, and decay constants for 
7T, p, K and K* mesons using our LFQM and compare 
with other theoretical model predictions. As expected, 
while the odd Gegenbauer moments for tt and p meson 



DAs are found to be zero due to isospin symmetry, the 
odd moments for K and K* meson DAs are nonzero due 
to the flavor SU(3)-symmetry breaking effect. We com- 
pute the decay constants for the transversely polarized 
p and K* mesons(/J and fx*) as well as the longitudi- 
nally polarized ones(/ p and fx* ) and compare with the 
light-cone sum rule(LCSR) calculations, in which the ra- 
tio of fy and fv is an important ingredient for the LCSR 
predictions of the B — > p and B — > K* transition form 
factors. We also confirm that our results for the d ecay 
constants follow an old SU(6) symmetry relation [39 |. 

The paper is organized as follows: In Sec. II, we briefly 
describe the formulation of our LFQM [3^, [3(| and the 
procedure of fixing the model parameters using the varia- 
tional principle for the QCD-motivated effective Hamilto- 
nian. The shape of the quark DA is then uniquely deter- 
mined in our model calculation. In Sec. Ill, the formulae 
for the quark DAs and decay constants of pseudoscalar 
and vector mesons are given in our LFQM. The Gegen- 
bauer and £(= x\ — x<i) moments are also given in this 
section. In Sec. IV, we present the numerical results for 
the decay constants, the quark DAs, the Gegenbauer and 
£ moments for (tt, K, p, K*) mesons and compare with 
other theoretical model predictions. Summary and con- 
clusions follow in Scc.V. The relations between £ and 
Gegenbauer moments are presented in Appendix A. 



II. MODEL DESCRIPTION 



by 



In our LFQM [35|, |36| , the meson wave function is given 



* J ^(x,k ± ,\\) = tf> R (x,k ± )n J J;{x,-k ± ), (4) 



where 0fl(x,k^) is the radial wave function and 
Tl^i* (x, kj_) is the spin-orbit wave function obtained by 
the interaction-independent Melosh transformation 
from the ordinary equal-time static spin-orbit wave func- 
tion assigned by the quantum numbers J PC . The meson 
wave function in Eq. is represented by the Lorentz- 
invariant variables Xi — pf /P + , k±i = p±i — XiP± and 
Ai, where P,Pi and are the meson momentum, the 
momenta and the helicities of the constituent quarks, re- 
spectively. 

The radial wave function </>/j(cc, kj_) of a ground state 
pseudoscalar meson (J^ = '") is given by 



(f> R (x,k±) 



1 



7T 3 / 2 /3 3 



1/2 



exp(-fc 2 /2/3 2 ), (5) 



where k 2 = k^ + k 2 and the gaussian parameter (3 is 
related with the size of the meson. Here, the longitudinal 
component k z of the three momentum is given by k z — 
(xi — \)Mq + (mf, — m 2 )/2Mo with the invariant mass 



M 2 



k 2 



k 2 



ml 



Xl 



X-2 



(6) 
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where x\ — x and X2 — 1 — x. The covariant form of the 
spin-orbit wave functions for pseudoscalar( J PC = h ) 
and vector(l ) mesons are given by 



/V AA 



^f2[Ml-{m 1 -m 2 )' i Yl^ 



TZ 1J - 3 = -- 



u(pi, A) 



Ah 



■(P1-P2) 



Mo+mi+ma 



«(P2, A) 



\/2[Af 2 - (mi -m 2 ) 2 ]V2 



(7) 



The polarization vectors & — (e + ,e used in this 

analysis are given by 



e*(±l) 



0,— e x (±)-P±,e±(±l) 



e±(±l) = T 



V2 : 



^(0) 



1 

Mo 



M 2 



(8) 



Note that V A} /^/^ = 
our wave function is given by 



1. The normalization of 



A A 



I d 3 k\^(x,k^XX)\ 2 



J\x J d 2 k ± (^)|<Mz,k 



x)| 2 = l, (9) 



where the Jacobian of the variable transformation 
{x, k±} -> k = (k_L,fc z ) is given by 



TABLE I: The constituent quark masses m q (in GeV) and the 
gaussian parameters /3 q q(m GeV) for the linear[HO] potential 
obtained from the variational principle. q=u and d. 



Pin 



0.22[0.25] 0.45[0.48] 0.3659[0.3194] 0.3886[0.3419] 



where Vconf = br(br 2 ) for the linear(HO) potential and 
(S q -Sq) — l/4(— 3/4) for the vector (pseudoscalar) meson. 

We then take (j>ii(x, kj_) as our trial function to mini- 
mize the central Hamiltonian via 



d(V\[H + V ]\y) 

d(3 



0. 



(13) 



From the above constraint, only 4 parameters are inde- 
pendent among the light-quark masses and the potential 
parameters {m q , j3 qq -, a,b, K)(q — u,d). In order to de- 
termine these four parameters from the two experimen- 
tal values of p and tt masses, we take the string tension 
b = 0.18 GeV 2 and the constituent u and d quark masses 
m u = m d = 0.22(0.25) GeV for the linear (HO) poten- 
tial, which are rather well known from other quark model 
analyses commensurate with Regge phenomenology [42j . 
More detailed procedure of determining the model pa- 
rameters of light quark sector(w(d) and s) can be found 
in [H, . Our model parameters for the light quark sec- 
tor obtained by the variational principle are summarized 
in Table [J 



III. QUARK DISTRIBUTION AMPLITUDES 
AND DECAY CONSTANTS 



dk z 



Mn 



dx 4x 1X2 



(mi — m 2 y 



Ml 



(10) 



The effect of the Jacobi factor has been analyzed in 
Ref. 

The key idea in our LFQM [3f| H(| for mesons is to 
treat the radial wave function 4>b.{x, kj_) as a trial func- 
tion for the variational principle to the QCD-motivated 
Hamiltonian saturating the Fock state expansion by the 
constituent quark and antiquark. The QCD-motivated 
effective Hamiltonian for a description of the meson mass 
spectra is given by [42j 



H qq =H + V qq = + fc 2 + ^/m 2 - + k' 2 + V qq . (11) 

In our LFQM [3^, we use the two interaction po- 
tential Vqq for the pseudoscalar and vector mesons: (1) 
Coulomb plus harmonic oscillator(HO), and (2) Coulomb 
plus linear confining potentials. In addition, the hyper- 
fine interaction, which is essential to distinguish vector 
from pseudoscalar mesons, is included for both cases, viz., 



4k 2S ■ S- 

V qq = V + Vh yp = a + Vconf - ^- + V -V 2 V C o 

6r 3mqm,q 



(12) 



The quark DA of a hadron in our LFQM can be 
obtained from the hadronic wave function by integrat- 
ing out the transverse momenta of the quarks in the 
hadron(see Eq. J2])), 



(x,fi) = 



|kx|<M d 2 k 



V167T 3 



r)k 

^*(x,kj.,AA). (14) 



For K and K* meson cases, we assign the momentum 
fractions x for s-quark and (1 — x) for the light u{d)- 
quark. The quark DA describes probability amplitudes 
to find the hadron in a state with minimum number of 
Fock constituents and small tranverse-momentum sepa- 
ration defined by an ultraviolet (UV) cutoff /1 > lGcV. 
The dependence on the scale /i is then given by the 
QCD evolution equation[lJ and can be calculated per- 
turbatively. However, the DAs at a certain low scale can 
be obtained by the necessary nonperturbative input from 
LFQM. Moreover, the presence of the damping Gaussian 
factor in our LFQM allows us to perform the integral up 
to infinity without loss of accuracy. The quark DAs for 
pseudoscalar (P) and vector(V) mesons are constrained 
by 



(j> P (y){x,n)dx 



fp(v) 
2V6 ' 



(15) 
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where the decay constant is defined as 

(0\q^J 5 q\P) = ifpP", (16) 
for a pseudoscalar meson and 

(0\q^q\V(P,X)) = fvMve^X), 
(0\qa^q\V(P,X)) = i/£[e"(A)P„ - e"(A)P M ], (17) 

for a vector meson with longitudinal(A = 0) and 
transverse(A = ±1) polarizations, respectively. The con- 
straint in Eq. (fT5"]) must be independent of cut-off [i up 
to corrections of order A 2 /u 2 , where A is some typical 
hadronic scale(< 1 GeV) [1(. For the non-perturbative 
valence wave function given by Eq. ([5]), we take fj, ~ 1 
GeV as an optimal scale for our LFQM. 

The explicit form of a pseudoscalar decay constant is 
given by 



where $> aa (x) = 6x(l — x) is the asymptotic DA and the 
coefficients a n (p) are Gegenbauer moments p| [l?], EH- 
The Gegenbauer moments with n > describe how much 
the DAs deviate from the asymptotic one. The zeroth 
Gegenbauer moment is fixed by the decay constant [l[, 
e.g. for the pion: 



a = 6 / [dx]4> v {xu n) 
'o 



V6 



(24) 



h 



tlx 



[d 2 k ± ] 



A 



^JA? 



0H(a;,kj 



(18) 



where A = (l—x)mi+xm2- The decay constants, fy and 
fy, for longitudinally and transversely polarized vector 
mesons, respectively, are given by 



fv= dx [d 2 kj 



ki 



f{- 



dx 



[d 2 k ± 



VA? 
(/>R(x,k±) 



A- 



2k 2 



Mq + mi + m 2 

(19) 



yM 2 



ki 



A 



The pion decay constant /£ xp ' 



M + mi + m 2 

(20) 

131 MeV is meaured 



from 7r — » pv and the p meson decay constant /p Xp ' — 215 

MeV is measured from p — ► e + e~ with the longitudinal ; 
polarization. While the constant jy are known from ex- ues for the ratio of h and fv 



where f„ ~ 131 MeV. In addition to the Gegenbauer 
moments, we can also define the expectation value of the 
longitudinal momentum, so-called ^-moments: 

<c> = J der*($ = J dxcHx), (25) 

where $(x) = 2$(2x — 1) normalized by = 1. In 
Appendix A, the relations between (£ n ) and a„(/j) are 
explicitly given up to n = 6. 



IV. NUMERICAL RESULTS 

In our numerical calculations, we use two sets of the 
model parameters for the linear and harmonic oscilla- 
tor(HO) confining potentials given in Table |TJ 

We show in Table HT1 our predictions for the decay con- 
stants of (tt, K, p, K*) mesons and compare with other 
theoretical model predictions [H, [3(| as well as data [H[ . 
As one can see from Table UH our results for the de- 
cay constants of (tt, K, p, K*) obtained from both lin- 
ear and HO potential models ( especially, those from HO 
potential) are compatable with the data UHl- Our val- 



i-e. fl/f p = 0.76[0.. 



periment, the constant fy are not that easily accessible 
in experiment and hence can be estimated only theoret- 
ically. 

The average value of the transverse momentum is given 

by 



(21) 



(K)qq= / «|ki||<Mx,kj 



Numerically, we have confirmed that (k^)^ = (3qq. 
This is a nonperturbative measure of the transverse size 
in the mesonic valence state. 

We may also redefine the quark DA as $>pry\(x) = 
(2>/6 / fp(y))<j)(x) so that 



<&p(v){x)dx = 1- 



(22) 



The quark DA <!>(x) evolved in the leading order (LO) 
of a s (p) is usually expanded in Gegenbauer polynomials 

n 3/2 



71=1 



a n (p)C 3 J 2 (2x-l) 



,(23) 



and fx./fK* = 0.82[0.86] obtained from the linear[HO] 
potential, are quite comparable with the recent QCD 
sum rule result^ fj/f p = (0.78 ± 0.08) and f^,/f K * = 
(0.78 ± 0.07) We also find that our results for the 

decay constants agree surprisingly well with an old SU(6) 
symmetry relation [HI], fJ( K .) = (U(K) + / P (/f*))/2 via 

the sum rule 0j (if ,) = (^(K)+^p(jf*))/ 2 > where ^ L(T) ( x ) 
is the longitudinally(transversely) polarized vector meson 
DA. 

We show in Fig. [T] the normalized quark DAs $p(a;) 
for 7r and K mesons obtained from linear (solid line) and 
HO (dashed line) potentials. For the pion DA, we also 
compare our results with the asymptotic result $ as (a;) = 
6x(l - x) (dotted line) as well as the A dS/CFT pre- 
diction $ AdS/CFT^) = tt\/x(1 — x)/8(double-dot- 
dashed line). For the pion case, our quark DAs obtained 
from both model parameters are somewhat broader than 
the asymptotic one. We also note from the normalized 
pion DAs that the suppression of the end point [x — > 
and 1) region has the following order in DAs, $tio{x) > 
$™(x) > <J>Linear(a:) > $Ads/CFT (x) ■ As discussed in 
Ref. [23|], there exists correlation between the shape of 
nonperturbative quark DA and the amount of low/high 
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TABLE II: Decay constants(in MeV) for the linear[HO] po- 
tential models compared with other models and data. 



/m Linear [HO] SRf30] Latticed] 



Exp. [45] 



U 

fK 
fp 
fp 
/*• 
fh 



130(131] 
161(155] 
246(215] 
188(173] 
256(223] 
210(191] 



205(9) 
160(10) 
217(5) 
170(10) 



126.6(6.4) 130.70(10) (36) 
152.0(6.1) 159.80(1.4) (44) 
239.4(7.3) 220(2) (a) , 209(4) (6) 



255.5(6.5) 



217(5) 



(c) 



(a) Exp. value for T(p° -> e+e~). 
O Exp. value for F(r — * pv T ). 
(c) Exp. value for V(r -> K*v T ). 

Q 2 contributions to the pion form factor. As the end- 
point region for the quark DA is more suppressed, the 
soft (hard) contribution to the pion form factor gets sup- 
pressed(enhanced). This finding is rather similar to the 
previous findings from the Sudakov suppression of the 
soft contribution [24, l25l . 1261 |27j. For the kaon case, the 
quark DA is asymmetric due to the flavor SU(3) symme- 
try breaking effect. The peak points of quark DAs for 
two potential models are moved slightly to the right of 
x = 0.5 point indicating that s-quark carries more lon- 
gitudinal momentum fraction than the light u(d)-quark. 

In the LO QCD the information of the leading- 
twist pion DA can be extracted from the pion-photon 
transition form factor F 7ri (Q 2 ) as follows: 



Q 2 F^(Q 2 ) 



V2U 



twist — 2 



\x^ M \. (26) 



o 6x(l — x) 



is 



The experimental value obtained in CLEO [2l| 
Q 2 F^(Q 2 ) = (16.7 ± 2.5 ± 0.4) x 10~ 2 GcV at Q 2 = 8 
GeV 2 , which goes to y/2f„ ^0.185 GeV in the asymp- 
totic Q 2 — > oo limit. With our leading twist pion DA 
shown in Fig. [TJ we obtain Q 2 F^(Q 2 ) = 0.202(0.181] 
GeV for the linear [HO] potential. 

For comparison between the leading twist and next- 
to-leading twist contributions to Q 2 F 7T ~ / {Q 2 ), we show in 
Fig. [Hour previous LFQM [31 prediction for Q 2 F^{Q 2 ) 
compared with the data [2l|, [44(. The thick solid and 
thick dashed lines represent our linear and HO potential 
model predictions including the higher twist effects(i.e. 
kj^ and the constituent mass m = m u = md) obtained 
from 



F 



NLO 



(Q 2 



(el 



-e 2 )^l 

7T3/2 J{i 

Mk 2 ) (1 



dx / d k 



dx 



x)m 



xjm 2 + k 2 ^ k' 



2 ' 



(27) 



where N c is the color factor and k'j_ = k^ + (1 — 
x)q±. The thin solid and thin dashed lines repre- 
sent our leading twist contribution(see Eq. (|26p) from 
^Linear i x ) and $ H o( a; ): respectively. We also com- 
pare our results with the leading twist contributions 



T 



T 



Linear potential(m=0.22 GeV, [5=0.3659 GeV) 

HO potential(m=0.25 GeV,p=0.3194 GeV) 

— WftM-Ix(1-x)|" 2 




0.2 0.4 0.6 



Linear potential 

HO potential 




FIG. 1: Normalized DAs $(x) for n and K mesons obtained 
from linear (solid line) and HO (dashed line) potential models 
compared with asymptotic result (dotted line) as well as the 
AdS/CFT prediction [H (double-dot-dashed line). 



from the asymptotic $ as (x) (dot-dashed line) and the 
AdS/CFT ^Ads/CFT^) (double-dot-dashed line). One 
should note that the AdS/CFT prediction^ \J x{\ — x)) 
increases the usual PQCD prediction^ x(l — x)) by 
4/3. Our higher twist results for both potential mod- 
els are not only very similar to each other but also 
in good agreement with the experimental data up to 
Q 2 ~ 10 GeV 2 region. At large Q 2 region, our higher 
twist prediction for the linear [HO] potential approaches 
Q 2 F^ LO (Q 2 ) = 0.194(0.180] compared to the leading 
twist result 0.202(0.181]. While the higher twist ef- 
fect on Q 2 F 7ri (Q 2 ) is large for the low and interme- 
diate Q 2 {< 10 GeV 2 ) region, its effect becomes very 
small for large Q 2 region compared to the leading twist 
contribution. Incidentally, it has been found that the 
leading Fock-state contribution to F^ 1 {Q 2 ) fail to re- 
produce the Q 2 = value corresponding to the axial 
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0.25 




0.2 










0.15 




(fit 




0.1 






— F^with^Cx) - 

F ^y LO With ^Linear (X> 


0.05 






— F^ 1 " with ^(x) 







, , 1 


" " F ky L ° with *Ads/CFr< x > 
♦ CLEO 
□ CELLO 



Q 2 [GeV 2 ] 

FIG. 2: The leading twist F^°(Q 2 ) and the next-to-leading 
twist F^y LO (Q 2 ) contributions to 7r — 7 transition form factor. 
Data are taken from Refs. [2ll. Ii3]. 



anomaly [13, [H[, i.e. it gives only a half of what is needed 
to get the correct ir° — > 77 rate [4{|. However, as shown 
in Refs. [13, HH H2] , the leading Fock-state contribution 
to F 7T1 (Q 2 ) has been enhanced by replacing the leading 
Fock-state wave function to an 'effective' valence quark 
wave function that is normalized to one. By taking the 
'effective' pion wave function with the asymptotic-like 
DAs, the authors in [HclIEHIEIJ found an agreement with 
the experimental data. Our LFQM prediction [35| also 
uses the same approach as Refs. [Ho|, HH, HI], i- e - tnc 
leading Fock-state 'effective' wave function that is nor- 
malized to one. The reason why our model is so suc- 
cesful for F 7r7 transition form factor is because the Q 2 
dependence(~ 1/Q 2 ) is due to the off-shell quark prop- 
agator in the one-loop diagram and there is no angular 
condition [53| associated with the pseudoscalar meson. 

In Tables HTT1 and ITVl we list the calculated (£ n ) and 
Gegenbauer moments a n (fi) for the pion (Tabic [Hi]) and 
the kaonfTable llVp DAs obtained from the linear[HO] po- 
tential models at the scale /i ~ 1 GeV. We also present 
the comparison with other model estimates at the scale 
of f < \i < 3 GeV. While the odd Gegenbauer mo- 
ments of the 7T meson DA become zero due to isospin 
symmetry, the odd moments for the kaon are nonzero 
due to a flavor-SU(3) violation effect of 0(m s — m u (d))- 
For the pion case, our result for the second Gegenbauer 
moment, aj = O.I2[0.05] obtained from linear[HO] po- 
tential is quite comparable with other theoretical model 
predictions given in Table A fair average is, how- 
ever, al = 0.17 ± 0.15 with still large errors [ll|. The 
LCSR based CLEO-data analysis @, i, 

on the tran- 
sition form factor suggests a negative value for aj, 
which is consistent with the result aJ(lGeV 2 ) > —0.07 
obtained in Ref. j|. Our result a\ = -0.003[-0.03] 
obtained from the linear [HO] potential also prefers a 
negative value consistent with the recent LCSR based 



CLEO-data analysis [ESQ- For the kaon case, the 
first moment is proportional to the difference between 
the longitudinal momenta of the strange and nonstrangc 
quark in the two-particle Fock component of the kaon, i.e. 
a i = (5/3)(x s — Xu) (x 3 — x,x B = 1 — x). The knowl- 
edge of af is important for predicting SU(3)-violation 
effects within any QCD approach that employs the quark 
DAs of mesons. In our model calculation, we obtain 
a positive value of af = 0.09[0.13] for the linear[HO] 
potential. Our results for af are quite consistent with 
those obtained from other estimates such as the previ- 
ous LFQM [23 (af = 0.08), the chiral-quark model [3 
(af = 0.096) and the QCD sum-rules [H, H|(af = 
0.05 ± 0.02). The positive sign of af can be understood 
intuitively since the heavier strange quark(antiquark) 
carries a larger longitudinal momentum fraction than the 
lighter nonstrange antiquark(quark). It is interesting to 
note that af for the HO potential model is greater than 
that for the linear one although the constituent mass dif- 
ference m s — m u = 0.23 GeV is the same for both models. 
This difference is attributed to the fact that the strange 
quark mass for HO potential model(m s = 0.48 GeV) is 
larger than for the linear one(m s = 0.45 GeV) and leads 
to more asymmetric shape for the HO potential. For 
the second Gegenbauer moment af , however, the linear 
potential model gives positive value(0.03), while the HO 
one gives negative value(-0.03). We note that the QCD 
sum-rules [H, [29| and lattice calculation [Tl[ give positive 
values while the chiral quark model [l5| gives a negative 
value. 

We show in Fig. [3] the normalized quark DAs for 
7r (upper panel) and K (lower panel) mesons obtained 
from the linear potential model(exact solution) and com- 
pare with those from the truncated Gegenbauer polyno- 
mials up to n = 6 (approximate solution). For the pion 
case, the truncation up to n — 4(dashed line) seems more 
close to our exact solution(solid line) than the truncation 
up to n = 6 (dot-dashed line) although the end-point be- 
havior of n = 6 case is more close to the exact solution 
than n < 4 case. Since the result up to n — 2 truncation 
is not much different from that up to n = 4 truncation, 
we do not show the result for n = 2 case in the figure. For 
the kaon case, while both truncations up to n — 4 (dashed 
line) and n = 6 (dot-dashed line) show good agreement 
with the exact solution(solid line), the truncation up to 
n = 2(dotted line) deviates a lot from the exact solution. 
Thus, it seems not sufficient to truncate the Gegenbauer 
polynomials only up to n = 2 for the kaon case. In both 
cases of tt and K mesons, our model calculation shows 
that the truncation of the Gegenbauer polynomials up to 
n = 4 seems to give a reasonable approximation to the 
exact solution. 

We show in Fig. 2] the normalized quark DAs Q p (x) 
for the longitudinally(solid line) and transversely(dashed 
line) polarized p meson obtained from the linear (upper 
panel) and HO (lower panel) potential models and com- 
pare with the asymptotic result(dotted line). For both 
potential models, the quark DA $^ with longitudinal 
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FIG. 3: Normalized DAs for 7r(upper panel) and A"(lower 
panel) mesons obtained from linear potential model compared 
with those obtained from the truncated Gegenbauer polyno- 
mials up to n = 6. 



polarization is not much different from the asymptotic 
result and the quark DA $ T p with transverse polarization 
is somewhat broader than both and <& as . Overall, 
the quark DAs for the p meson are closer to the asymp- 
totic result than those for the pion case. Although the 
overall shapes of our $^ and <1>J are not much different 
from the asymptotic result, the end-point behaviors of 
our model calculation exhibiting the concave shape are 
different from the asymptotic result of the convex shape. 
We also should note that our DAs for the p meson satisfy 
the SU(6) symmetry relation [H, $J = ($„. + $£)/2. 

In Table [V] we list the calculated (£") and Gegen- 
bauer moments a n (p) for the p meson DAs obtained 
from the linear [HO] potential models at p ~ 1 GeV and 
compare with other model estimates. The odd Gegen- 
bauer moments of both longitudinally and transversely 
polarized p meson DAs become zero due to isospin sym- 



metry. Our £ moments obtained from both linear and 
HO potential models are similar to the asymptotic re- 
sults and quite consistent with the previous LFQM [2(| 
and QCD sum-rules [3(j. However, slight differences of 
£ moments among different model predictions turn out 
to be quite sensitive in terms of Gegenbauer moments. 
For instance, the second Gegenbauer moment a^ L for the 
longitudinally polarized p meson obtained from the lin- 
ear potential model gives positive value(0.02), while the 
same from the HO potential model gives negative value(- 
0.02). This may be still comparable with the previous 
LFQM [20|_giving negative value(-0.03) and the QCD 
sum-rules [30| giving positive value^.O^Q^y). Our re- 
sults for the higher Gegenbauer moments a p n (n > 4) ob- 
tained from both linear and HO potential models show 
negative values regardless of polarization states. These 
negative values for the higher Gegenbauer moments are 
related with the concave shapes of quark DAs, <I>^ (x) and 
at the end-point region. 

We show in Fig. [5] the normalized quark DAs <&k * {%) 
for longitudinally(solid line) and tranversely(dashed line) 
polarized K* meson obtained from the linear(upper 
panel) and HO (lower panel) potential models. As in the 
case of the p meson, the shape of for the transversely 
polarized K* meson near the central(.T = 1/2) region is 
somewhat broader than that of for the longitudi- 
nally polarized K* meson in both linear and HO models. 
Also, the peak points for both <I>^.» and are shifted 
to the right of x = 1/2 point due to the SU(3) flavor sym- 
metry breaking as in the case of K meson. Our quark 
DAs for the K* meson satisfy also the SU(6) symmetry 
relation jjHj, § T K » = (<f> K + $%.)/2. 

In Tables El and EH we list the calculated (f l ) 
and Gegenbauer moments a n (p) for the longitudi- 
nally(Table IVI[) and transversely(Table IVII[) polarized 
K* meson DAs obtained from the linear [HO] potential 
models at p, ~ 1 GeV and compare with the available 
QCD sum-rule results (U HJ. While the odd Gegen- 
bauer moments of the p meson DAs become zero due to 
the isospin symmetry, the odd moments for the K* me- 
son DAs are nonzero because the SU(3) flavor symmetry 
is broken. Our values of the first Gegenbauer moments, 
af£ = 0.11 [0.14] and a% = 0.10[0.14] for the linear[HO] 
potential model are in a good agreement with the QCD 
sum-rule results [§cj af L = a^ T = 0.10±0.07 at the scale 
p = 1 GeV. Note that the positive af^ refers to K* con- 
taining an s quark but the sign will change for K* with 
an s quark. Our predictions for the even powers of a^nh 
and af 2 \-i'(n > 1) give negative values while the LCSR. re- 
sults [30] give positive values. However, the recent QCD 
sum-rule result af T = 0.02 ± 0.02 [3l| at the scale p = 1 
GeV is consistent with our value a^ T = 0.04[0.03] ob- 
tained from the linear[HO] potential model. Also, the 
(£ 2n ) moments are not much different between the two 
models. We show in Fig. [6] the normalized quark DAs 
for transversely polarized p(upper panel) and if* (lower 
panel) mesons obtained from the linear potential model 
and compare with those from the truncated Gegenbauer 
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FIG. 4: Normalized DAs for the longitudinally (solid line) and 
transversely(dashed line) polarized p meson obtained from lin- 
ear(upper panel) and HO (lower panel) potential models com- 
pared with asymptotic result(dotted line). 



polynomials up to n = 6. For the p meson case(upper 
panel), since the truncation up to n — 2 does not much 
differ from that up to n = 4 (dashed line), we do not 
show the result for n = 2 case in the figure. Both trunca- 
tions up to n — 2 and n — 4 are quite close to our exact 
solution(solid line). The truncation up to n = 6 (dot- 
dashed line) having a deep at x = 1/2 point shows a 
slight deviation from the exact solution. For the K* me- 
son case(lower panel), both truncations up to n = 4 and 
n = 6 show good agreement with the exact solution, 
while the truncation up to n = 2 (dotted line) deviates 
a lot from the exact solution. For both p and K* meson 
cases, the truncation of the Gegenbauer polynomials up 
to n = 4 seems to give an overall reasonable approxima- 
tion to the exact solution. 



Linear potential model 

- * L K -(x) 

- — * T K <x) 




x 



FIG. 5: Normalized DAs for longitudinally (solid line) and 
tranversely(dashed line) polarized K* meson obtained from 
linear(upper panel) and HO(lower panel) potential models, 
respectively. 



V. SUMMARY AND DISCUSSION 

In this work, we investigated the quark DAs, the 
Gegenbauer moments, and decay constants for it, p, K 
and K* mesons using the LFQM constrained by the vari- 
ational principle for the QCD-motivated effective Hamil- 
tonian. Our model parameters obtained from the varia- 
tional principle uniquely determine the above nonpertur- 
bative quantities. 

Our predictions for the quark DAs for it and p mesons 
show somewhat broader shapes than the asymptotic one. 
The odd Gegenbauer moments for 7r and p meson DAs 
become zero due to isospin symmetry. Our predic- 
tions for aj and a\ are consistent with the recent light- 
cone sum-rule based CLEO data analysis for the pion- 
photon transition form factor. Interestingly, we also find 
that our leading twist result for the pion-photon tran- 



9 



4> [1 + up to n=4| 
4> [1 + up to u=6| 
4>, M (x)(Linear potential) 




tf> (x)| 1 + up to 0=2] 
3> (x)[l + up to n=4| 
3> (x)[l + up to n=6| 
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FIG. 6: Normalized DAs for the transversely polarized 
p(upper panel) and K* (lower panel) mesons obtained from 
linear potential model compared with those obtained from 
Gegenbauer polynomials up to n — 6. 



sition form factor, Q 2 F^{Q 2 ) = 0.202[0.181] GeV ob- 
tained from the linear [HO] potential model, is reduced 
to Q 2 F^ LO (Q 2 ) = 0.194[0.180] if we include the higher 
twist effects such as the transverse momentum and the 
constituent mass. Our result is quite compatible with 
the CLEO data, Q 2 F 7ri (Q 2 ) = (16.7 ± 2.5 ± 0.4) x 1CT 2 
GeV at Q 2 = 8 GeV 2 [2l|. The quark DAs for K 
and K* show asymmetric forms due to the flavor SU(3)- 
symmetry breaking effect. This leads to the nonzero val- 
ues of the odd Gegenbauer moments. In our model cal- 
culations of the quark DAs for (tt, K, p, K*) mesons, the 
truncation of the Gegenbauer polynomials up to n = 4 
seems to give a reasonable approximation to the exact 
solution. 

Our predictions for the decay constants for 7T, K, longi- 
tudinally polarized p and K* mesons are in a good agree- 
ment with the data. We also obtain the decay constants 



for the transversely polarized p and K* mesons (/J and 
/#») using our LFQM. Our predicted values of fy/fv 
averaged between the linear and HO potential cases 
are (fJ/f P U = 0.78 and tf%./f K *)«, = 0.84. They 
are consistent with the recent QCD sum rule results, 
fj/fp = (0.78 ±0.08) and /£.//*. = (0.78 ±0.07) ^. 
Moreover, our results for the decay constants are in a 
good agreement with the SU(6) symmetry relation [39l |. 

= (A-(-ff) + f P (K*))/2 via the sum rule <% x .) = 
(^•n-(if) + 0p(_R"*))/2- Further investigations to utilize our 
LFQM are underway. 
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APPENDIX A: RELATION BETWEEN £ AND 
GEGENBAUER MOMENTS 



The ^-moments (£") defined by Eq. ([25|) can be related 
to the Gegenbauer moments a n (p) in Eq.J22|). The rela- 
tions up to n — 6 are given by 



if) 
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12 

35 
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35 
3 
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35 21 J ' 
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Also, the first six Gegenbauer polynomials in Eq. 



(Al) 
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as follows: 
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(2ie 4 - uf + i), 

£(33£ 4 -30£ 2 + 5), 
(3003£ 6 - 3465£ 4 + 945£ 2 - 35). (A2) 
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TABLE III: The £ and Gegenbauer moments cCCaO f° r the pion DAs obtained from the linear[HO] potential models compared 
with other model estimates. The numbers in the parentheses stand for the scales of the corresponding works. 



Models 


r> 


(O 


(O 


al 


a\ 


< 


Linear [HO] 


0.24[0.22] 


0.11[0.09] 0.07[0.05] 


0.12[0.05] 


-0.003[-0.03] 


-0.02[-0.03] 


Asymptotic WF 


0.20 


0.09 


0.05 











AdS/CFT [22J 


0.25 


0.125 


0.078 


0.146 


0.057 


0.031 


[3] (2.4 GeV) 


0.35 


0.21 




0.44 


0.25 


- 


[4] (1.0 GeV) 


0.24 


0.11 




0.115 


-0.015 


- 


[5] (1.0 GeV) 


0.28 


0.13 


_ 


0.23 


-0.05 


- 


[8] (1.0 GeV) 


0.27 


0.12 




0.20 


-0.14 


- 


[14] (1.0 GeV) 


0.22 


0.10 




0.046 


0.007 


- 


[15](1.0 GeV) 


0.21 


0.09 


0.05 


0.029 


-0.046 


-0.019 


[16] (1.0 GeV) 


0.22 


0.09 




0.05 


-0.04 




[11] (2.0 GeV) 


0.269(39) 






0.201(114) 






[6] (2.4 GeV) 


0.24±0.01 






0.12±0.03 






[8] (2.4 GeV) 


0.25 


0.11 




0.14 


-0.08 




[10] (2.4 GeV) 


0.23 


0.11 




0.08 


0.02 




[21] (2.4 GeV) 


0.27 


0.11 




0.19 


-0.14 




[46] (2.4 GeV) 


0.25 


0.12 




0.16 


0.02 





TABLE IV: The £ and Gegenbauer moments a% (fi) for the kaon DAs obtained from the linear[HO] potential models compared 
with other model estimates. The numbers in the parentheses stand for the scales of the corresponding works. 



Models 


W) 


r> 


<n 


(O 


(e) 




Linear [HO] 


0.06[0.08] 


0.21[0.19] 


0.03[0.04] 


0.09[0.08] 


0.02[0.03] 


0.05[0.04] 


[15] (1.0 GeV) 


0.057 


0.182 


0.023 


0.070 


0.012 


0.0345 


[20] (1.0 GeV) 
[29] (2.0 GeV) 


0.046 
0.03±0.01 


0.20 
0.26±0.04 


0.025 


0.08 


0.015 


0.04 


Models 








af 


ag 


a 6 K 


Linear [HO] 
[15] (1.0 GeV) 


0.09[0.13] 
0.096 


0.03[-0.03] 
-0.051 


0.06[0.04] 
-0.008 


-0.02[-0.03] 0.007[-0.01] 
-0.040 -0.002 


-0.01[-0.01] 
-0.0097 


[20] (1.0 GeV) 


0.08 





0.03 


-0.06 


-0.14 


-0.03 



[28](1.0 GeV) 0.05±0.02 0.27±S^ 
[11] (2.0 GeV) 0.0453(9)(29) 0.175(18)(47) 
2!) 12.0 GeV) 0.05±0.02 0.17+0.10 



TABLE V: The £ and Gegenbauer moments a£(/i) for the p meson DAs obtained from the linear[HO] potential models compared 
with other model estimates. The numbers in the parentheses stand for the scales of the corresponding works. 



Models 














Linear [HO] 
Asymptotic WF 


0.21[0.19] 
0.20 


0.09[0.08] 
0.09 


0.05[0.04] 
0.05 


0.22[0.20] 
0.20 


0.10[0.09] 
0.09 


0.06[0.04] 
0.05 


[20] (1.0 GeV) 
[30] (1.0 GeV) 


0.19 
n 90+0.03 

u - zo -0.02 


0.07 

11+ 003 

u - 11 -0.02 


0.036 


0.2 

n 90+0.03 

u -^°-0.02 


0.082 
11+ 003 

u - 11 -0.02 


0.042 


Models 


a 2L 


a iL 








"6T 


Linear [HO] 
[20] (1.0 GeV) 


0.02[-0.02] 
-0.03 


-0.01[-0.03] 
-0.09 


-0.02[-0.02] 0.06[0.007] 
0.7 


-0.01[-0.03] 
-0.04 


-0.02[-0.02] 
-0.04 


[30] (1.0 GeV) 


0.09i°;£° 


0.03±0.02 




0.091°,;" 


0.03±0.02 





TABLE VI: The £ and Gegenbauer moments a„ (/J-) for the longitudinally polarized K* meson DAs obtained from the linear[HO] 
potential models compared with other model estimates. The numbers in the parentheses stand for the scales of the corresponding 
works. 



Models (^) L 






(C)l 






Linear [HO] 0.07 [0.09] 
[30](1.0 GeV) 0.06±0.04 


0.19[0.18] 

n 99+O.O3 
u - zI ^-0.02 


0.03[0.04] 


0.08[0.07] 
0.10±0.02 


0.02[0.02] 


0.04[0.03] 


Models af^ 


K" 
a 2L 


K~ 

a 3 L 


K~ 
0-4L 


K~ 
a 5L 


K' 


Linear[HO] 0.11 [0.14] 


-0.03[-0.07] 


0.03[0.02] 


-0.02[-0.03] 


0[-0.01] 


-0.01[-0.01] 


[30](1.0 GeV) 0.10±0.07 


n7 +009 

u - u ' -0.07 




0.02±0.01 
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TABLE VII: The £ and Gegenbauer moments a„ (fj,) for the transversely polarized K* meson DAs obtained from the linear [HO] 
potential models and compared with other model estimates. The numbers in the parentheses stand for the scales of the 



Models {£ L ) T 




(e 3 >T 






(C)t 


Linear [HO] 0.06 [0.08] 
[30] (1.0 GeV) 0.06±0.04 
[31](1.0 GeV) 0.06±0.04 


0.20[0.18] 

n 99+O.O3 
U.^_ .02 


0.03[0.04] 
0.03±0.02 


0.08[0.07] 
0.10±0.02 


0.02[0.02] 


0.04[0.04] 


Models afx 










a& 


Linear[HO] 0.10[0.14] 
[30](1.0 GeV) 0.10±0.07 
[31](1.0 GeV) 0.10±0.07 


-0.008[-0.05] 
ny+o.09 


0.04[0.03] 
0.02±0.02 


-0.02[-0.03] 
0.02±0.01 


0.003[-0.01] 


-0.01[-0.01] 



